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Abstract 

We consider the generalized continuous-time Lyapunov equation: 



H ■ A*XB + B*XA^-Q, 



where QisanNx N Hermitian positive definite matrix and A, B are arbitrary Nx N matrices. Un- 
der some conditions, using the coupled fixed point theorem of Bhaskar and Lakshmikantham, we 
establish the existence and uniqueness of Hermitian positive definite solution for such equation. 
Moreover, we provide an iteration method to find convergent sequences which converge to the 
solution if one exists. 

Keywords. Lyapunov equation, positive definite solution, iterative algorithm, coupled fixed point, 
(— I ' ordered metric space. 
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1 Introduction 

Consider the generalized continuous-time algebraic Lyapunov equation (GCALE) 
(N 

OO ; A*XB + B*XA=-Q (1) 



cn ! v«th given matrices Q,A,B and unknown matrix X. Such equations play an important role in stability 

O ■ theory [6, 18], optimal control problems [14, 17] and balanced model reduction [16]. In [19], it is 

proved that Eq. (1) has a unique Hermitian, positive definite solution X for every Hermitian positive 
definite matrix Q if and only if all eigenvalues of the pencil XA - B are finite and lie in the open left 
half-plane. When A = I (the identity matrix), Eq. (1) becomes the standard Lyapunov equation 

X' 

C^ ; XB + B*X = -Q. (2) 

The classical numerical methods to solve Eq. (2) are the Bartels- Stewart method [1], the Hammarling 
method [10] and the Hessenberg-Schur method [9]. An extension of these methods to solve Eq. (1) 
with the assumption A is nonsingular, is given in [5, 7, 8, 9, 19]. Other approaches to solve Eq. (1) are 
the sign function method [4, 13, 15], the ADI method [3, 11, 20]. 

In this paper, we consider Eq. (1), where Q is an N x N Hermitian positive definite matrix and 
A,B are arbitrary N x N matrices. Using Bhaskar-Lakshmikantham coupled fixed point theorem [2], 
we provide a sufficient condition that assures the existence and uniqueness of a Hermitian positive 
definite solution to Eq. (1). IVIoreover, we present an algorithm to solve this equation. Numerical 
experiments are given to illustrate our theoretical result. 
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2 Notations and preliminaries 

We shall use the following notations: M{N) denotes the set of all N x N matrices, HiN) c MiN) the 
set of all Nx N Hermitian matrices and P{N) c H(_N) is the set of allNx N positive definite matrices. 
Instead of X e P(_N), we will also write X > 0. Furthermore, X > means that X is positive semidefi- 
nite. As a different notation for X - F > (X - F > 0) we will use X > F (X > F). If X, F e HiN) such 
that X < F, then [X, F] will be the set of all Z e H[N] satisfying X < Z < F. If X, F e HiN) such that 
X < F, then ]X, F[ will be the set of all Z e HiN) satisfying X<Z<Y. 

We denote by || • || the spectral norm, i.e., \\A\\ = \/^^iA*A), where A"*" (A* A) is the largest eigenvalue 
of A* A. The N x N identity matrix will be written as /. 

It turns out that it is convenient here to use the metric induced by the trace norm || • || i . Recall that 
this norm is given by ||A||i = L/=i SjiA), where sjiA), j = l,...,n are the singular values of A. In fact, 
we shall use a slight modification of this norm. For QeP(Ar), we define ||A||i_q = ||Q^'^AQ^'^||i. 

For any U e MiN), we denote by Sp{[/) the spectrum of U, that is, the set of its eigenvalues. 

The following lemmas will be useful later. 

Lemma 2.1 (See [21]) LetA> andB >ObeNxN matrices, thenO< triAB) < \\A\\ ■ triB). 

Lemma 2.2 (See [12]) LetAe Hin) satisfying-I <A<I, then \\A\\ < 1. 

Definition 2.1 Let (A, <) be partially ordered set. We say that a mapping F : A x A ^ A has the mixed 
monotone property if for all X, Y,J,K e A, 

X<J,Y>K^> Fix, Y) < FiJ,BO. 

The proof of our main result is based on the following fixed point theorems. 

Tiieorem 2. 1 (Biiaskar and Laksimiikantiiam [2] ) Let iA.,<) bea partially ordered set and dbea met- 
ric on X such that (A, d) is a complete metric space. LetF : A x A ^ A foe a continuous mapping having 
the mixed monotone property on A. Assume that there exists a 5 e [0, 1) with 

5 

diFix,y),Fiu,v)) < -[dix,u) + diy,v)] 

for allx>u and y<v. Suppose also that 

(i) there exist Xq, ya e A such thatxo < Fixo,yo) andyo > F(yo,xo); 

(ii) every pair of elements has either a lower bound or an upper bound, that is, for every ix,y) e Ax A, 
there exists azeA such thatx < z and y<z. 

Then, there exists a unique x e A such thatx = F(x,x). Moreover, the sequences {x„} and {y„} defined 
byxn+i = Fix„,yn) andy„+i = F(y„,x„) converge tox, with the following estimate 

_ 6" 

max{dixn,x),diyn,x)] < -max{d(xo,Xi),(i(yo,yi)}, for all n. 

1-0 

Tiieorem 2.2 (Sciiauder Fixed point tiieorem) LetS beanonempty, compact, convex subset of a normed 
vector space. Every continuous function f :S^ S mapping S into itself has a fixed point. 
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3 Main result 

Our main result is tiie following. 

Theorem 3.1 Suppose that there exists (Q, M) e PiN) x PiN) such that 

(a) 2{UQU*+BQB*)<Q; 

(b) 2{VQV* +WQW*)<Q; 

(c) U*MU + B*MB<M-iV*MV+W*MW]; 

(d) Qe]U*MU + B*MB,M-iV*MV+W*MW)l, 

where 

A-B+I A+B+I 
U= — , V= — and W = B-I. 

Then, 

(i) Eq. (1) has one and only one solution X e P{N) . 

(ii) X elQ - [U* MU + B* MB),Q + [V* MV + W* MW)]. 

(iii) Let (X„) and (F„) the sequences defined by Xo = 0, Fq = M, and 

Xn+i = Q + [V*X„V + W*XnW)-{U*YnU+B*Y„B) 
F„+i = Q+{V*Y„V+W*YnW]-{U*XnU + B*XnB) 



We have 



(3) 



lim \\X„ - X\\ = lim || Y„ - X\\ = 0, (4) 

n— ►oo n—'oo 



and the error estimation is given by 

max{||X„-X|M|y„-X||}<Y3^max{||Xi-Xoll,||Fi-Fol|}, (5) 

for all n, where < 5 < 1 . 

Proof. It is easy to show that Eq. (1) is equivalent to 

X=Q + {V*XV+W*XW)-{U*XU + B*XB). (6) 

Consider the continuous mapping F : H{N) x H{N) — ' H{N) defined by 

F{X,Y) = Q + {V*XV+W*XW)-iU*YU + B*YB), forall X,YeHiN). (7) 

Clearly, Eq. (6) is equivalent to 

X=F[X,X). (8) 

Let X, YJ,Ke H[N] such that X<JandY>K. Then, 

F{X,Y) = Q + [V*XV + W*XW)-{U*YU + B*YB) 
< Q+[V*JV + W*JW)-[U*KU + B*KB) 
= F{J,K). 

This implies that F is a mixed monotone mapping. 
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Let X, Y,J,Ke H{N] such that X>JandY<K.We have 

\\F[X, Y)-F{J,K)\\^Q = ti{Q^'^[F[X, Y)-F{J,K))Q^'^) 

= tx[Q^''^{V*{X-})V + W*{X-nW+U*{K-Y)U + B*{K-Y)B)Q^''^] 

= lx[Q^''^iV\X-nV + W*{X-T)W)Q^''^) + tx[Q^''^{U*{K-Y)U + B*{K-Y)B)Q^'^). 

On the other hand, using Lemma 2.1, we have 

lx[Q^''^{V\X-})V+W*iX-J)W)Q^''^]=tx[VQV*{X-}) + WQW*{X-J)) 

= tx[VQV* {X- J)Q^''^Q-^'^ + WQW* {X- J)Q^''^Q-^''^] 

= tx[Q-^'^VQV*{X-J)Q^'^ + Q-^''^WQW*{X-})Q^''^) 

= lx[Q-^'^VQV*Q-^'^Q^'^[X- J)Q^'^ + Q-^'^WQW*Q-^'^Q^'^[X- J)Q^'^) 

<\Q-^''^VQV*Q-^'^^Q-^''^WQW*Q-^''^\tx[Q^''^{X-J)Q^'^) 

= \\Q~^'^[VQV* + WQW*)Q-^'^X.x[Q^''^{X- J)Q^'^] 

= \\Q~'^'^[VQV* + WQW*]Q-^'^\\ ||X-/||i_Q. 

Thus, we have 

tx{Q^'^[V*[X-J)V + W*[X-J)W)Q^'^]<\\Q~'^'^iVQV* +WQW*)Q~'^'^\\\\X-J\\^Q. (9) 
Similarly, we have 

tx[Q^'^{U*{K- Y)U + B*{K- F)B)Q1'2] < \\Q-'^'^[UQU* +BQB*)Q-^'^\\ \\K- Y\\^q. (10) 

Now, using (9) and (10), we get 

\\F[X,Y]-F[J,K]\\^q<-[\\X-J\\^q + \\K-Y\\^q], 

where 

5 = 2xxiax{\\Q-^'^{VQV* + WQW*)Q-^'^\\ , \\Q-^'^{UQU* + BQB*)Q-^'^\\} . 

Now, from Lemma 2.2, (a) and (b), we have 5 e (0, 1). 

Taking Xq = and Yq = A^> from (c) and (d), it follows that 

Xo<F(Xo,Fo) and Yo>FiYo,Xo). (11) 

Since all the hypotheses of Theorem 2.1 are satisfied, we deduce that there exists a unique X e 
HiN) solution to Eq. (8). This implies that Eq. (1) has a unique solution X e H{N). 

Now, to establish (i), we need to prove that X e P{N). The Schauder fixed point theorem will be 
useful in this step. Define the mapping G : [F(0, M),F[M, 0)] -^ H[N) by 

G{X)=F{X,X), forall Xe [F(0,M),F(M,0)]. 

Note that from (11) andthemixedmonotonepropertyof F, wehaveF(0,M) <F(M,0). 

We shall prove that G([F(0,M), F(M,0)]) c [F(0,M],F(M,0)]. LetXe [F(0,M],F(M,0)], that is, 

F(0,M)<X<F(M,0). 

Using the mixed monotone property of F, we get 

F(F(0,M],F(M,0)) <F(X,X) = G(X) < F(F(M,0),F(0,M)). 

On the other hand, from (11), we have 

0<F(0,M) and M>F(M,0). 
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Again, using the mixed monotone property of F, we get 

F{F{M,0),F{0,M))<F{M,0) and F(F(0,M),F(M,0)) > F(0,M). 

Then, 

F(0,M)<G(X)<F(M,0). 

Thuswe proved that G([F(0,M],F(M,0)]) c [i7(0,M],F(M,0)]. 

Now, G maps the compact convex set [F(0, M),F{M, 0)] into itself. Since G is continuous, it follows 
from Schauder's fixed point theorem (see Theorem 2.2) that G has at least one fixed point in this set. 
However, fixed points of G are solutions of (1), and we proved already that (1) has a unique Hermitian 
solution. Thus this solution must be in the set [F{0,M),F(M,0)], that is, 

Xe [Q- {U* MU + B* MB),Q + {V* MV +W* MW)](^ P{N). 
Thus, we proved (i) and (ii). The proof of (iii) follows immediately from Theorem 2.1. D 

Now, we present some consequences following from Theorem 3.1, when A,B are Hermitian ma- 
trices. 

Corollary 3.1 Suppose that 

(a) ABeHiN); 

(b) 2iUQU + BQB)<Q; 

(c) 2{VQV+WQW)<Q, 

where 

A-B+I A+B+I 
U= — , V= — and W = B-I. 

Vz Vz 

Then, 

(i) Eq. (1) has one and only one solution X e P{N) . 

(ii) X elQ -ZiUQU + BQB),Q + 2[VQV + WQW)]. 

(iii) Let{Xn) and {¥„) the sequences defined by Xq = 0, Fq = 2Q, and 

Xn+i = Q+iVX„V + WX„W)-{UYnU + BY„B) 
F„+i = Q+ iVYnV +WY„W) - {UXnU + BX„B) ' 



We have 



lim ||X„-X||= lim ||F„-X||=0, 

«— ►oo n—oo 



and the error estimation is given by 

max{||X„ -XII, II Y„ - X\\} < A_max{||Xi - XqII, || Y, - Fol|} , 
for all n, where < 5 < 1. 
Proof. It follows from Theorem 3.1 by taking Q = Q and M = 2Q. D 

Corollary 3.2 Suppose that 
(a) A,BeH{N); 
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(b) 2{U^ + B^)<I, 2{V^ + W^)<I; 

(c) U^+B^<Q<I-iV^ + W^), 

where 

A-B+I A+B+I 

U= — , V^ — and W = B-I. 

\/2 VZ 

Then, 

(i) Eq. (1) has one and only one solution X e PiN) . 

(ii) XelQ-iU^ + B^),Q + [V^ + W^)]. 
(iii) Let (_Xn) and (F„) the sequences defined by Xq = 0, Fq = /, and 

Xn+i = Q+iVXnV+WX„W)-[UYnU + BY„B) 
Yn+i = Q+iVY„V+WY„W)-iUX„U + BX„B) ' 



We have 



lim ||X„-X||= lim ||F„-X||=0, 



and the error estimation is given by 

max{||X„-l||,||F„-l||}<^max{||Xi-XolM|Fi-Foll}, 
for all n, where < 5 < 1. 
Proof. It follows from Theorem 3.1 by taking Q = M = I. D 

4 Numerical experiments 

In this section, we present some numerical experiments to check the convergence of the proposed 
algorithm (3). We take Xq = and Fq = M e P{N). For each iteration i, we consider the residual errors 

Ei [X) = II Xi - (Q + [V* Xi V+W* Xi W) - iU* Xi U + B*XiB))\\, 

EiiY) = II Yi - {Q + {V* YiV+W* YiW) - ([/* YiU + B* F,5))|| 

and 

Ei = manEi{X),Ei{Y)}. 

All programs are written in MATLAB version 7.1. 
Example 4.1 We consider Eq. (1) with 

' 2 0.02 0.05 \ 
Q= 0.02 2 0.02 , A=\ 
, 0.05 0.02 2 j 

In this case, we have A, Be //(3), 

Sp(Q-2([/Qf/+5Q5)) ={0.3345,0.3379, 0.3880}, Sp(Q-2(yQy+VKQVK)) = {o.4532,0.4573,0.4612}, 

which imply that conditions (a)-(c) of Corollary 3.1 are satisfied. 

Considering the iterative method (3) with Xq = and Yq = 2Q, after 100 iterations one gets an 
approximation to the positive definite solution X and it is 

I 1.9495 0.0288 0.0142 

X ~ Xioo = Fioo = 0.0288 1.9496 0.0288 

, 0.0142 0.0288 1.9495 

a«rf£ioo = 1.9215x10"^^. 



-0.95 


0.001 


0.001 \ 




( 0.54 


-0.002 


-0.002 


0.001 


-0.95 


0.001 


, B = 


-0.002 


0.54 


-0.002 


0.001 


0.001 


-0.95 j 




I -0.002 


-0.002 


0.54 
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Example 4.2 We consider Eq. (1) with 



Q = 



( 0.4 0.01 0.02 0.03 0.04 ^ 

0.01 0.4 0.01 0.02 0.03 

0.02 0.01 0.4 0.01 0.02 

0.03 0.02 0.01 0.4 0.01 

i^ 0.04 0.03 0.02 0.01 0.4 j 



A = 



'' -0.95 0.001 0.001 0.001 
0.001 -0.95 0.001 0.001 
0.001 0.001 -0.95 0.001 



0.001 "^ 

0.001 

0.001 



0.001 0.001 0.001 -0.95 0.001 
0.001 0.001 0.001 0.001 -0.95 



and 



0.44 


-0.02 


-0.02 


-0.02 


-0.02 ^ 


-0.02 


0.44 


-0.02 


-0.02 


-0.02 


-0.02 


-0.02 


0.44 


-0.02 


-0.02 


-0.02 


-0.02 


-0.02 


0.44 


-0.02 


-0.02 


-0.02 


-0.02 


-0.02 


0.44 



B = 



In this case, we have A, Be H{5), 

Spil - 2(f/2 + B^)] = |o.4078, 0.4078, 0.4078, 0.4078, 0.64716| , 

Spil -2{V^ + W^)] = |o.0094, 0. 1577, 0.1577, 0. 1577, 0.1577| , 

Spio - [U^ + B^)] = |o.0516, 0.0882, 0.0963, 0.0984, 0.3049| , 

and 

Spil-{V^ + W^) -Q] = |o.0231, 0.1844, 0.1865, 0.1949,0.23121, 

which imply that conditions (a)-(c) of Corollary 3.2 are satisfied. 

Considering the iterative method (3) with Xq = and Yq = I, after 82 iterations one gets an approx- 
imation to the positive definite solution X and it is 



X~Xi 



and Es2 = 7.0549 X 10 



82 ■ 



-16 



t82 



f 0.4895 0.0429 0.0541 0.0658 0.0781 ^ 

0.0429 0.4878 0.0418 0.0535 0.0658 

0.0541 0.0418 0.4873 0.0418 0.0541 

0.0658 0.0535 0.0418 0.4878 0.0429 

0.0781 0.0658 0.0541 0.0429 0.4895 
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